The effect of isoscalar-isovector coupling in infinite nuclear matter 
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Working on the framework of Relativistic Mean Field theory, we exposed the effect of nonlinear isoscalar- 
isovector coupling on G2 parameter set on the density dependence of nuclear symmetry energy in infinite nuclear 
matter. The observables like symmetric energy and few related coefficients are studied systematically. We 
presented the results of stiff symmetry energy at sub-saturation densities and a soft variation at normal densities. 
Correlation between the symmetric energy and the isoscalar-isovector coupling parameter fully demonstrated 
for wide range of density. The work further extended to the octet system and showed the effect of coupling over 
the equation of state. 

PACS numbers: 21.65.Ef, 21.10.Gv, 21.65.+f, 26.60.+C 



I. INTRODUCTION 

The best possible and well defined theoretical laboratory 
to studied for many body system is infinite nuclear matter at 
certain conditions. To review the status of microscopic stud- 
ies of nuclear matter and neutron rich matter, that already 
reaches to the destination by mean-field models (relativis- 
tic and non-relativistic) tjj-|6|] and although by some other- 
microscopic, methods such as Brueckner-Hartree-Fock (BHF) 
j7Ki3l and Dirac-Brueckner-Hartree-Fock (DBHF) ITlUll. 
The isospin and density dependence of the nuclear symmetry 
energy, E sym , is the current interest because of its implica- 
tions not only in nuclear physics but also in astrophysics. The 
softening in equation of state (EOS) of nuclear matter likely 
lead to an exciting problem in astro-nuclear physics from few 
decades to till date. Mean while the novel phenomena like 
formation of superheavy nuclei in astrophysical system is the 
front learning about the island of stability for that regions. 

The prediction of compressional modulus {K ) in symmet- 
ric nuclear matter for nonrelativistic and relativistic models 
differ by about 25 in magnitude. In other word, the nonrela- 
tivistic models predict Kq ~ 220-235 MeV Jj-lg], relativistic 
models argue for K ~ 250-280 MeV GUI- The density de- 
pendence of the symmetry energy E sym is not exposed fully at 
pres ent, which may be responsible for the above discrepancy 
01611 . Note that the E sym is the difference between the en- 
ergy of pure neutron matter and normal nuclear matter, which 
is a tremendous approximation tfH [13. Ell • The large value 
of Kq shows appropriately stiff symmetric energy i.e. rises 
rapidly with baryon density 016111711 . This make a passing ref- 
erence between Kq and E sym , which is a function of density 
and isospin component of scalar and vector mesons. Thus, 
it is enlightening to explore the effects of E sym in the RMF 
model, which was limited to a narrow range from the analysis 
of skin data, on the composition and structure of hot proton 
neutron and cold neutron star matter that hold a large density 
range JHH. 

In this work, we introduce an extra term on the top of the 
Gl or G2 forces to the Lagrangian, which is the combined 
effect of scalar and vector field with isoscalar- isovector cou- 
pling constant (A„). The newly added term affect to the sym- 



metric energy and the slope parameter L sym . The problem 
is the stiffness of the E sym with respect to baryonic density 
followed a softer path with respect to density along with this 
extra added coupling. Without this additional coupling con- 
stants one may not be overcome the hindrance in the nuclear 
matter. The most important prospectus of the coupling con- 
stant is that the various other properties and nuclear matter 
observables are not seriously compromising in their magni- 
tude. 



II. THE RELATIVISTIC MEAN FIELD (RMF) THEORY 

The complete description of EOS can be implemented by 
using quantum chromodyanmics (QCD), which is the funda- 
mental theory of the strong interaction. The concept of effec- 
tive field theory (EFT) is applicable to low ener gy Q CD lf24ll 
and is called as quantum hadrodynamics (QHD) B25M29I1 . The 
description of normal nuclear matter properties II25L 13011 and 
finite nuclei JUrH^l are successfully explained by QHD. The 
exchange of mesons like iso-scalar scalar a, iso-scalar vector 
ui, iso-vector vector p and the pseudoscalar it are responsible 
for nucleon interaction and increases the degree of freedom. 
The nucleons are considered as Dirac spinors moving in clas- 
sical meson fields. In mean-field approximation, the contribu- 
tion of the 7T meson is zero, because of its pseudospin nature. 
The extension of the standard relativistic mean-field (RMF) 
theory with addition of non-linear scalar-vector and vector- 
vector self interaction is called as E-RJVIF proposed by Furn- 
stahl, Serot, and Tang If3~5ll36ll (see also Miller and Serot If37ll . 
The coupling constants are self-consistent with the underlying 
symmetries of QCD by non-renormalization. The concept of 
naturalness can focus by expanding the non-linear Lagrangian 
and organized it in increasing power of field and their deriva- 
tive according to their dimensional analysis 11381 13911 and trun- 
cated at a given level of accuracy 04Oh42I1 . Here we have taken 
the fourth order of expanded field and in addition to this we 
have taken the special interaction between the u and p with 
a new coupling constant A„ to get good result in the study 
of E sym and other co-efficients. The modified Lagrangian is 
given as JHQ: 
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The co-variant derivative D p is defined as 



-Dp = dp + ig^p + ig P hT a p 



a n a 



(2) 



where R a and f2 p „ are field tensors and defined as follow 



(3) 



where, = \/(fc 2 + m b* 2 ) is the effective energy, k B is the 
Fermi momentum, rn* B = niB—ga^ stands for effective mass 
(which solved self-consistently), I3 is the third component of 
isospin projection and B stands for baryon. The energy and 
pressure density depends on the effective mass m* B of the sys- 
tem, so first need to solve this self-consistent Eqn. (4). As the 
self-consistent field obtained, the expression for energy and 
pressure density are given as; 
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and 



(4) 



Here ms is the mass of nucleon, and m a , and m p are 
the masses of corresponding mesons. The equation of motion 
derived from the Lagrangian using mean-field approximation 
self-consistently. The field equation for different meson fields 
are given by 
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(ii) 



where Pi and are lepton's pressure and energy density, re- 
spectively (and here we have taken electron only). Note that 
the following quantities are redefined as, Hq = g u u>o and 
R = g p po3- At zero temperature, the pressure of the system 
obtained from energy density and its first derivative, which is 
thermodynamically consistent in mean-field (MF) approxima- 
tion. 
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III. CALCULATIONS AND DISCUSSION 



The scalar (p s ) and vector (p) densities are given as follow: 
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We have solved the mean field equations self-consistently 
and parametrized the coupling constant by renormalization. 
We have made attention to G2 l35l l36ll parametrization sets 
of the E-RMF model to reproduce the fundamental proper- 
ties of nuclear matter. In this parameter set £ coupling con- 
stant is determined by self-interactions of w-meson, which is 
responsible for the softening the EOS at high density and re- 
produced the maximum mass of neutron star. For softening 
the symmetric energy of the system at nuclear matter density, 
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FIG. 1: The EOS or energy per nucleon and pressure density as a 
function of density from RMF compare with the other predictions. 



100 



75 

> 

i 

^50 



■— ■ RMF (G2) 

□— o GSkll 

O—O SQMC650 

lSkT2 
<J— <1 Ska35s20 
H — ^ SKRA 
O— O Skxs20 




FIG. 2: The results of symmetric energy E sym is a function of den- 
sity p from RMF for different values of A„ compare with others pre- 
dictions. 



we have added an extra coupling between iso-scalar and iso- 
vector field with coupling constant A„. This made make pos- 
sible to soften the E sym and reproduced saturation properties 
for the force parameter at A„ ~ 0.00. In other word, A„ is sim- 
ply a self interaction constant between isoscalar and isovector 
field, which does not change any other nuclear matter proper- 
ties like energy density, pressure density and so many consid- 
erably. 



B. Symmetry Energy 

The symmetric energy E syrn , which is important in infinite 
nuclear matter and finite nuclei because of isospin dependent 
in the interaction. The expression of this key observable is 
defined as 
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A. Energy and Pressure density 



where the effective energy E* F — \/{k F + to* 2 ), kp is the 
Fermi momentum and the effective mass m* = to — gs<fio- 
The effective mass of the p-meson modified because of non- 
linear coupling (p — lo) interaction and is given by 



Infinite nuclear matter is important for the investigation of 
physical quantities relevant to heavy nuclei and compact ob- 
jects like neutron star. At saturation density po, the binding 
energy per particle e as a function of density is well estab- 
lished quantities consistent to thermodynamic relation at cer- 
tain temperature. These translational invariance of the sys- 
tem facilitates theoretical calculations of these observables. 
The prediction of E-RMF (G2+A„) for EOS of nuclear matter 
(NM) (solid black line) compared with DBHF fl45[] and other 
theoretical predictions l46ll are in Fig. 1 (a). From the fig- 
ure, it is clear that all the prediction well matches with respect 
to density. Fig. 1(b), shows the comparison of pressure den- 
sity with experimental data l47ll and theoretical predictions 
lf45L l48l l49h . The obtained results for energy density (see Fig. 
1) from Ref. ll2lll is exactly same to the present calculation. 
This implies that the EOS of SNM is not affected by the new 
coupling with A„. Being insensitive to this observable, the 
aim of this present work to pursue the systematic variation of 
E sym by employing the isoscalar-isovector coupling parame- 
ter A„, which discussed below. 
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where m p is the mass of the p-meson. It is noted that the non- 
linear isoscalar-isovector coupling (A„) modified the density 
dependent of E sy m without affecting the saturation properties 
of the SNM (as shown in Fig. 1 (a)). In this model the sym- 
metric nuclear matter saturates at a Fermi momentum of = 
1.00 /to 1 with E ~ 16.0 MeV and an compressibility of Kq 
= 215 MeV ll2lll . We can expand the density dependence of 
the Esym with the help of these 3-co-efficients as a function 



of y — (p — p ) /3po as follow: 
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FIG. 3: The slope parameter L aym along with K sym and Q aym for different value of Av, with G2 parameter set. 



K S ym = 3 po L sym . (16) 

V dp J P=P0 

and 

Q aym = 27pl(^p) . (17) 

Here L sym is called the slope parameter, which is nothing but 
the slope of E sym . The quantity K sym represents as the cur- 
vature of E sym with respect to density. In SNM, the value of 
K syrn at saturation density (po) entitled by Kq and called as 
symmetry incompressibility. Finally, Q sy m is the third deriva- 
tive of the symmetric energy, which called as skewness coef- 
ficient [46]. As shown in the Fig. 2, we plotted the E sym 
for RMF (G2+A t ,) compare with other theoretical predictions 
j46Tl . The softness of the E aym with the coupling A v clearly 
observed from the figure. Further, the rate of change of E sym 
i.e. the slope parameter L sym for different values of A„ is dis- 
played in Fig. 3. Both the figure shows similar behaviour with 
the coupling constant and with baryonic density. 

There are many references which trying to fix the con- 
straint on symmetric energy and L sym l46l |50T - I55ll . Finally 
we comes to the constraint on E sym and L sym at saturation 
density (po) and the obtained results are compared with with 
experimental data at HIC Q, PDR |H3 HI, IAS JH and 
FRDM [60] in Fig. 4 and Table 2. From the figure, it is clear 
that for A„=0.00, our observables are not within the exper- 
imental constraint region. For the values of A„=0.01-0.04, 



TABLE I: The results for slope co-efficient S sym , K sym for different 
value of A„ are in list. 



A„ 


Esym 


Lsym 


0.00 


36.40 


105.67 


0.01 


38.41 


111.69 


0.02 


38.87 


113.09 


0.03 


38.62 


112.32 


0.04 


38.02 


110.52 


0.05 


37.27 


108.27 


0.06 


36.46 


105.85 


0.07 


35.64 


103.42 


0.08 


34.86 


101.05 


0.09 


34.11 


98.79 


0.10 


33.39 


96.65 



these comes within limit and again go out with further in- 
creasing the value of A„ =0.05. Hence, from this behavior, 
we may suggest the limiting value of A„ is 0.01-0.04, which 
can improved the observables with G2 force parameter. 



IV. HADRON MATTER 

The density at the interior core of the neutron star is very 
high - 7-8/9 so other hadronic part IH l6lU6l also 
come to play role on the equation of state. In this section, 
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ates at high density. The octet EOS compares well with the 




FIG. 4: The Symmetric energy So with respect to L aym at saturation 
density compare with flow data. 
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we extended our calculations to the system with octet (n, p, 
A°,£ 0,±1 , The purpose of this work to show the effect 

of isoscalar-isovector interaction with coupling constant A„ 
on the equation of state of the nuclear system with full octet. 
It is worthy to mention that, many authors l2ll l4ll Ir54ll have 
been used Gl and G2 forces for the description of finite and 
infinite nuclear matter properties. Moreover, the non-linear 
isoscalar-isovector interaction in mean field was discussed in 
Refs. Ir35[l66ll . here we suppose to check this interaction with 
G2 along with A„, keeping all saturation constraints. We as- 
sume that all the hyperon in the octet system have the same 
coupling ratio with mesons. This is given by a constant fac- 
tor which is equal to y/2/3 It54ll . the density of the system ~ 
7-8 po and assuming that the stars are the composite system 
of baryons. The electron e _1 and muon \x are for maintai ning 
the charge neutrality and [}— equilibrium condition llrSil l67ll 
under the week interactions: 



Si -> S 2 + 1 + W, B 2 +l-^B 1 



n, 



where Si, S2, I, v, v are baryons, leptons, neutrino and antin- 
uetrino respectively. In case of octet system, we are dealing 
with neutron star (real system), which became unstable with 
small change in the system. The equation of state obtained 
for nuclear matter with nucleon and hadrons for G2 + A v is 
shown in Fig 5. 

We plotted pressure(/m -4 ) density with respect to energy 
density (/m~ 4 ) in Fig. 5. Here we can see that at high den- 
sity octet EOS becomes stiffer by increasing the A„ value but 
it behave different way at low density region. Equation of 
state for necleon case (EOS-NP) becomes softer with increas- 
ing A„ value over whole the density region. We compare our 
results with the empirical EOS obtained by Steiner et cd. Ir58ll . 
Along with the octet we also plotted nucleon EOS at different 
A„ values. Nucleon EOS well matches with the data of Steiner 
et al. 16811 at low and high density but the octet EOS is devi- 



FIG. 5: Equation of states(EOS) with the E-RMF model for octet and 
NP at different A„ value. 



experimental data only at intermediate density. An impressive 
observation is noted from this plot, the influence of A„ with 
G2 parameter set is not much because in G2 already taken 
the four power of iso-scalar meson so it is more effective then 
iso-scalar and iso-vector interaction. 



V. SUMMARY AND CONCLUSIONS 

After analyzing the result whatever got from this calcula- 
tion, we forcefully say that A„ i.e. the coupling of isoscalar 
and isovector interaction is very crucial to softening the equa- 
tion of state of nucleonics and octet system. One can eas- 
ily observed that the E sym curve with respect to p/po became 
softer with varying A„ but after A„ ~ 0.05 it became constant. 
Same result we got in L sym which is the slope of symmetric 
energy. For symmetric nuclear matter isoscalar and isovector 
interaction do not effect the E/A and pressure density curve 
with respect to p/po- We got identical curve as in Ref. I2TI1 
for Gl and G2 sets with experimental data 84711 . In case of 
hadron, the extra term containing A„ plays an important role 
in softening the equation of state. Hence, we may say that this 
interaction term is crucial for both cases (nucleon and hadron). 
Further, we extend our calculation to the neutron star proper- 
ties in near future. 
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